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CHAPTER  1 
INTRODUCTION 

1.1.  Preview 

This  study  considers  optimal  decision  strategies  with  regard  to 
capacity  expansion  in  a stochastic  demand  environment.  Chapter  1 is  an 
introduction  to  the  topic  at  hand  and  provides  the  preliminaries  upon 
which  later  model  construction  is  based.  As  Indicated  in  the  next  section, 
there  are  two  types  of  facilities  which  must  be  considered  in  a capacity 
expansion  model:  permanent  and  temporary.  With  regard  to  temporary 

facilities,  a further  classification  is  necessary  according  to  whether 
or  not  the  temporary  facilities  are  modular. 

Model  I,  presented  in  [7],  treats  the  case  where  temporary 
facilities  are  non-modular.  Model  II,  presented  herein,  treats  the  case 
where  temporary  facilities  are  modular.  For  both  cases,  it  is  shown  that 
the  determination  of  optimal  expansion  sizes  can  be  expressed  as  a single- 
variable,  nearly-unconstrained,  minimization  problem  of  a very  particular 
form.  The  solution  to  this  problem  is  treated  in  [8],  where  a revised 
model  is  also  introduced  to  demonstrate  a number  of  generalizations  that 
are  possible  in  both  Models  I and  II. 


1.2.  Permanent  and  Temporary  Facilities 

When  one  considers  the  growth  in  capacity  of  a singular  system, 
such  a3  a plant,  pipeline,  superhighway  or  school,  a common  phenomena 
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often  arises.  Specifically,  when  demand  first  exceeds  the  initial  capacity 
of  the  system,  a permanent  capacity  expansion  is  not  immediately  under- 
taken. Instead,  temporary  operating  measures  are  instigated  in  order  to 
satisfy  excess  demand  over  the  short-run.  In  some  cases,  the  temporary 
measure  taken  may  be  simply  to  backlog  excess  demand.  In  other  cases, 
where  backlogging  is  undesirable,  the  present  system  may  be  overloaded 
beyond  its  planned  capacity  in  order  to  accommodate  excess  demand.  In 
still  other  cases,  where  backlogging  is  undesirable  and  overloading  is 
infeasible,  temporary  capacity  may  be  rented  from  outside  sources.  In  order 
to  distinguish  between  the  normal  accommodation  of  demand  and  the  emergency 
accommodation  of  excess  demand  through  temporary  measures,  facilities 
will  be  categorized  here  as  either  permanent  or  temporary.  That  is, 
"temporary  facilities"  will  denote  the  measures  utilized  to  accommodate 
excess  demand  over  the  short-run  prior  to  the  implementation  of  a full- 
fledged  expansion  of  permanent  facilities. 

The  use  of  temporary  facilities  usually  involves  some  sort  of  cost 
penalty  over  that  of  permanent  facilities;  otherwise,  one  might  well  ask 
which  facilities  are  really  temporary.  Thus,  the  economic  advantages  of 
temporary  facility  usage  might,  on  the  surface,  seem  questionable.  There 
are,  however,  two  principal  factors  which  promote  the  desirability  of 
temporary  facility  usage.  Firstly,  there  is  the  problem  of  uncertainty: 
when  demand  first  exceeds  the  permanent  capacity,  is  the  excess  a short-term 
peak  or  a long  term  trend?  If  the  time  interval  of  excess  is  to  be  of 
small  duration,  then  an  expansion  of  permanent  facilities  can  be  disastrous, 
since  overhead  costs  on  unused  capacity  can  easily  lead  to  financial 
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insolvency.  Secondly,  there  is  the  problem  of  capitalization:  permanent 

capacity  expansion  costs  often  exhibit  significant  economies  of  scale, 
which  in  turn  dictate  substantial  expansion  sizes  requiring  large  capital 
expenditures.  Even  if  sustained  demand  growth  is  certain,  bankruptcy  is 
possible  over  the  short-run  if  sufficient  revenue  cannot  be  generated  from 
new  facilities  to  meet  capitalization  costs.  Furthermore,  capitalization 
itself  may  prove  difficult  until  excess  demand  is  of  sufficient  magnitude 
to  convince  investors  that  permanent  capacity  expansion  is  warranted. 

While  each  of  the  above  problems  can  in  itself  promote  the  use  of 
temporary  facilities,  the  two  taken  together  generate  an  even  stronger 
case  for  temporary  facility  usage  under  an  expected  discounted  cost 
criterion.  Figure  1. 1 provides  some  illustration.  The  demand  peaks 
at  times  t^  and  t 2 are  short-lived;  any  large  permanent  expansion 
at  these  times  will  result  in  a great  deal  of  wasted  capacity  over  a 
significant  portion  of  the  time  interval  shown.  Temporary  measures 
can  be  used  to  accommodate  these  excesses.  Similarly,  at  time  t^,  the 
long-term  prospects  for  demand  are  unclear  and  temporary  facilities  can 
be  used.  However,  by  time  t^,  the  excess  demand  has  grown  to  significant 

size  so  as  to  somehow  "warrant"  a permanent  expansion  (at  which  time,  the 
temporary  facilities  in  use  would  be  discontinued).  If  temporary  facilities 
are  not  used,  then  a permanent  expansion  must  be  undertaken  at  time  t^. 
Thus,  the  use  of  temporary  facilities  forestalls  the  permanent  expansion 
for  (t^-t^)  time  units.  If  costs  are  discounted,  then  the  discount 
savings  alone,  accruing  from  the  prolongment  of  the  permanent  expansion 
capital  outlay,  can  easily  exceed  any  cost  penalties  resulting  from  the 
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FIGURE  1.1:  Use  Temporary  Facilities  to  Forestall  Permanent  Expansion 


permanent  capacity.  When  k is  nonnegative,  the  permanent  facilities 
suffice  to  serve  all  demand.  However,  when  k is  negative,  excess 
demand  exists  and  temporary  facilities  must  be  used.  Let  K denote  the 


limit  on  temporary  facility  usage.  Whenever  k = -K  (that  is,  excess 
demand  equals  K)  and  a new  unit  of  demand  growth  occurs,  a permanent 
expansion  of  size  X+l  (X  > K)  is  required.  When  permanent  expansion 
occurs,  the  temporary  facilities  in  use  are  discontinued  and  the  spares 
level  k increases  to  X-K  > 0.  If  demand  consists  of  deterministic 


constant  growth,  then  the  use  of  this  type  of  recurrent  (X,K)  expansion 
policy  results  in  a saw-tooth  pattern  for  the  spares  level,  as  illustrated 
in  Figure  1.2. 


The  value  K serves  as  a "trigger",  forcing  a permanent  expansion 
whenever  a new  unit  of  demand  growth  occurs  while  k = -K.  Thus,  K+l 
represents  the  point  at  which  excess  demand  warrants  a permanent  expansion. 
In  cases  of  overloading,  the  value  K may  be  determined  based  on  physical 
constraints  governing  the  degree  of  overloading  that  can  safely  be  incurred. 
In  other  cases,  K may  be  determined  according  to  a judgment  concerning 
the  willingness  of  investors  to  participate  in  capitalization  and  the 
prospects  of  generating  sufficient  revenue  to  meet  capitalization  costs, 
given  that  excess  demand  has  reached  the  value  X.  In  such  cases  where 
K Is  predetermined,  the  parameter  of  interest  will  be  X*(K)+1,  the 
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FIGURE  1.2: 


Recurrent  (X,K)  Expansion  Policy  for  Linear  Demand 
Spares  Level  = Permanent  Capacity-Demand. 
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optimal  permanent  expansion  size  to  be  utilized  in  conjunction  with  the 
temporary  facilities  usage  limit  K.  In  other  cases,  a value  K*  will 
be  determined  based  on  a cost  minimization  over  all  feasible  operating 
policies  (X,K).  Given  a procedure  for  determining  X*(K),  this  last 
problem  reduces  to  that  of  a cost  minimization  over  all  feasible  operating 
policies  (X*(K),  K). 


1.3.  Modular  Temporary  Facilities 

Modular  temporary  facilities  can  be  viewed  as  a specific  case 
of  jobletting  where  temporary  facilities  can  be  rented  from  outside 
sources  in  a fixed  increment  size.*  An  instance  of  this  type 
la  the  use  of  mobile  "barracks"  facilities  (e.g.,  house  trailers)  as 
auxiliary  classrooms  to  alleviate  conjestion  in  overcrowded  schools. 

Another  instance  of  this  type  is  the  use  of  "pair-gain"  devices 
to  augment  the  capacity  of  over- subscribed  telephone  cables.  To  illustrate, 
consider  Figure  1.3.  In  order  to  simplify  the  illustration,  consider  only 
the  one-way  communication  from  a subscriber  telephone  to  the  switching 
equipment  at  the  subscriber's  local  telephone  office.  As  shown  in  Figure 
1.3a,  this  communication  is  normally  implemented  by  the  dedication  of  a 
single  wire-strand  from  the  subscriber  telephone  to  the  office.  This 
single  wire  is  typically  one  strand  of  a telephone  cable  between  the 
subscriber's  neighborhood  and  the  office.  Suppose  that  the  demand  for 

*Note : Modular  cost  behavior  can  also  be  found  in  cases  of  backlogging 
and  overloading;  however,  it  is  more  prevelant  in  jobletting  or 
rental  situations. 
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telephones  in  this  neighborhood  increases  to  a point  at  which  all  wire- 
strands  of  the  existing  cable  arc  used.  The  telephone  company,  being 
required  to  provide  service  to  all  soliciting  subscribers,  must  somehow 
increase  the  cable  capacity  in  order  to  satisfy  new  demands.  One  obvious 
possibility  is  the  immediate  placement  of  a new  cable  (i.e.,  an  immediate 
permanent  expansion).  Another  possibility  is  the  implementation  of 
multiple- party  service;  however,  this  alternative  is  usually  unexceptable 
to  consumers  and  is,  in  fact,  not  used  in  many  areas.  Another  alternative 
to  the  placement  of  an  additional  cable  is  the  utilization  of  pair-gain 
devices  on  the  existing  cable.  As  illustrated  in  Figure  1.3b,  a pair-gain 
device  consists  of  electronics  that  permit  the  transmission  of  a number 
of  simultaneous  conversations  over  a single  wire-strand  by  encoding  the 
conversations  at  one  end  and  decoding  the  conversations  at  the  other  end. 

As  shown  in  the  figure,  these  devices  are  typically  available  in  a fixed 
modular  size  allowing  the  transmission  of  L additional  conversations 
over  a single  wire- strand;  the  "pair-gain  ratio"  of  each  such  module  is 
said  to  be  "L  to  1". 

The  cost  characteristic  which  distinguishes  modular 
temporary  facilities  from  non-tnodular  temporary  facilities  is  the  existence 
of  instaneous  charges  over  and  above  the  normal  rental  costs.  Specifically, 
an  installation  charge  is  typically  incurred  whenever  an  additional 
module  is  engaged  and  similarly,  a removal  charge  is  incurred  whenever 
a module  is  returned  to  the  outside  supplier.  Thus,  in  the  case  of 
modular  temporary  facilities,  there  is  a sub-optimization  problem 
with  regard  to  how  modules  should  be  engaged  and  returned  in  order  to 
avoid  excessive  installation  and  removal  charges. 
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1.4. 


The  Poisson  Demand  Model 


In  order  to  account  for  the  promotion  of  temporary  facility  usage 
due  to  uncertainty^  demand  will  be  recognized  as  being  stochastic  in 
nature.  Specifical ly;  it  will  be  assumed  that  demand  increases  ("arrivals") 
are  characterized  by  a Poisson  process  [ 6 ] at  rate  A-^  > 0.  Similarly,  it 
will  be  assumed  that  demand  decreases  ("departures")  are  characterized  by 
a Poisson  process  at  rate  Xg  > 0.  Furthermore,  the  arrival  and  departure 
processes  will  be  presumed  to  be  independent  of  each  other  and  also 
independent  of  the  expansion  policy  chosen. 

Given  the  above  demand  charac ter ization,  define  an  "event"  as 
either  an  arrival  or  a departure.  Then  an  equivalent,  and  more  convenient, 
characterization  is  as  follows:  events  occur  according  to  a Poisson 

process  at  rate  A = A^  + A0;  the  probability  that  an  event  is  an  arrival 
equals  p = Aj/X;  and  the  probability  that  an  event  is  a departure  equals 
q = Ag/X  ( *1  = 1-p)  L b ].  That  is,  the  demand  process  treated  here  can 
be  viewed  as  a random  walk  with  exponential  inter-event  times.* 

It  should  be  noted  that  the  above  demand  characterization  is 
pointed  toward  systems  providing  a homogeneous  continuous  service  to  a 
fairly  stable  clientele,  such  as  a power  plant  or  a telephone  service  area. 

It  is  really  not  intended  to  model  actual  c lien  tele  movement  for  cases 
where  customers  are  rapidly  flowing  in  and  out,  such  as  a job  shop.  In 
these  cases,  it  is  suggested  that  "customer  demand"  be  viewed  as  the 
maximum  usage  level  required  of  facilities  during  short  intervals  of 
time  (e.g.,  weeks  or  months). 

★ 

Notice  that  the  departure  process  is  untruncated : even  when  the  demand 
level  is  zero,  departures  are  allowed.  The  viability  of  this 
assumption  is  demonstrated  in  Section  1.7  of  [7]. 
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It  should  also  be  noted  that  the  demand  process  Is  assumed  time- 
stationary. For  this  reason,  the  system  being  modelled  must  be  in  the 
midst  of  a relatively  stable  period  of  either  sustained  growth  or  sustained 
negative  growth. 

Finally,  notice  that  the  arrival  and  departure  rates  are  presumed 
constant  with  regard  to  the  demand  level.  From  a practical  point  of  view, 
this  is  convenient  since  only  two  estimates  --  a single  arrival  rate  and 
a single  departure  rate  --  are  necessary  to  implement  the  model.  However, 
in  most  realistic  situations,  it  appears  that  the  arrival  and  departure 
rates  will  vary  according  to  the  system  demand  level.  In  recognition  of 
this  fact,  the  revised  model  of  [ 8 ] permits  these  rates  to  vary. 

1.5 . Costs 

The  optimization  criterion  used  here  will  be  the  minimization  of 
all  future  expected  discounted  costs;  r will  denote  the  continuous  dis- 
count rate,  0 < r < 1.  Three  types  of  costs  will  be  allowed:  permanent 

expansion  costs,  permanent  facility  operating  charges  and  temporary 
facility  charges.  All  costs  are  assumed  to  be  time-stationary. 

Permanent  expansion  costs  will  be  denoted  by  g(*),  where  g(X) 
is  the  cost  of  a permanent  expansion  of  size  X+l.  Notice  that  g(*) 
is  presumed  to  be  a function  only  of  the  expansion  size,  independent  of 
the  level  of  existing  permanent  capacity  and  the  value  of  the  temporary 
facilities  usage  limit,  K.  This  restriction  is  relaxed  in  the  revised 
model  of  [ 8 ] . 
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It  1 9 important  to  note  that  g(0)  represents  the  cost  of  a unit 
expansion.  Thus,  g(0)  is  comprised  of  the  fixed  (i.e.,  "setup")  expan- 
sion cost  plus  the  marginal  cost  for  the  first  unit  of  increased  capacity. 
Hence,  any  presumption  regarding  a specific  functional  form  (e.g.,  con- 
vexity) for  g over  [(),<»)  will  not  preclude  the  existence  of  a fixed 
expansion  cost. 

With  regard  to  other  costs,  the  temporary  facility  charges  are  of 
principal  interest  in  ►his  study.  In  order  to  provide  sufficient  generality 
with  regard  to  these  costs,  the  permanent  facility  operating  charges 
will  be  presumed  to  be  proportional  to  the  amount  of  permanent  capacity  used . 
The  following  theorem  indicates  the  modelling  simplifications  that  result 
fr cm  this  assumption. 

Theorem  1.1.  Suppose  that  permanent  facility  operating  charges  are  pro- 
portional to  the  amount  of  permanent  capacity  used,  at  rate  y per  unit- 
time. Then  an  equivalent  cost  model  is  given  as  follows: 

( i)  When  permanent  facilities  satisfy  all  demand,  no  operating  charges 
are  incurred. 

(ii)  When  temporary  facilities  are  necessary  (i.e.,  k < 0)  marginal 

costs  equal  to  the  temporary  facility  charges,  less  an  adjustment 
y|k|  = -yk  > 0 per  unit-time,  are  incurred. 

Proof . See  Theorem  1.1  of  | 7 ) . 
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The  presumed  independence  of  the  permanent  operating  costs  from 
the  permanent  facility  capacity  may  not  be  that  restrictive,  since  pro- 
portional fixed  operating  charges  which  depend  solely  on  the  permanent 
capacity  level  can  be  included  in  the  expansion  cost  function  g.  To 
illustrate,  suppose  that  the  actual  permanent  operating  costs  include 
a proportional  fixed  charge  of  per  unit- time  for  each  unit  of  per- 
manent capacity  available.  Let  YQ  denote  the  initial  permanent  capacity. 
The  expected  discounted  cost  resulting  from  the  ijf  charge  on  the  initial 

oo  rt  ^ 

capacity  is  then  / YQ\)f  e*  dt  = \|rr  Y^j  this  cost  is  independent  of 
the  expansion  policy  chosen.  The  expected  discounted  costs  resulting 
from  the  \(r  charge  on  future  capacity  expansion  increments  can  then  be 
accounted  for  by  defining  g as 

g(X)  = g (X)  + / (X-t-1) \|r  e"rt  dt 

0 

= ge(X)  + Vr_1(x+l)  , 

where  gg  represents  the  actual  expansion  function  alone. 

Temporary  facility  costs  will  be  represented  by  a per-module 
rental  charge  n per  unit-time,  installation  charge  c,  and  removal 
charge  d;  notice  that  the  later  two  costs  are  instantaneous  charges. 

All  restrictions  (except  time-stationarity)  on  the  facility 
charges,  both  permanent  and  temporary,  are  also  relaxed  in  the  revised 
model  of  ( 8 ] . 
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1.6.  Remarks  and  Notation 


Related  results  are  discussed  at  length  in  Section  1.6  of  [7], 
Pertinent  related  models  include  those  of  Manne  |5J,  Freidenfelds  [2J, 
and  Koontz  [4].  (7]  also  demonstrates  a near-analogy  between  the  model 

at  hand  and  single-item  inventory  models  with  both  returns  and  backlogs. 

As  introduced  previously  in  this  chapter,  £( c)  will  denote  the 
demand  at  time  t > 0,  X will  denote  the  event  rate  with  arrival  prob- 
ability p and  departure  probability  q,  k will  denote  the  spares 
level,  K will  denote  the  temporary  facilities  usage  limit,  X+l  will 
denote  the  expansion  size  and  g will  denote  the  permanent  expansion 
cost  function,  where  g(X)  represents  the  cost  for  an  expansion  of 
size  X+l. 

IR111  will  denote  the  space  of  all  real  vectors  having  m com- 


ponents, m > 1.  Both  row  and  column  vectors  will  be  utilized;  however, 

0)  T 

a vector  is  a column  vector  unless  otherwise  noted.  Given  V 6 IS  , V 

mxn 

will  denote  the  transposed  vector.  R will  denote  the  space  of  all 
real  matrices  having  m rows  and  n columns.  Given  A G ERm  n,  A^ 


will  denote  row  i of  A and  similarly, 
j of  A (Af  € IRn  and  A G IR®). 

1 • • J 


A . will  denote  column 


CHAPTER  2 


MODEL  II:  THE  MODULAR  CASE 

2. 1.  Introduction 

This  chapter  characterizes  the  untruncated  Poisson  model  under  two 
basic  assumptions: 

Assumption  2.1.  Permanent  facility  operating  charges  are  proportional 
to  the  amount  of  permanent  capacity  used,  at  rate  per  unit-time. 

Assumption  2.2.  Modular  temporary  facilities  are  available  to  satisfy 
excess  demand.  Each  module  provides  L units  of  capacity.  The  installa- 
tion of  a module  incurs  an  instantaneous  connect  charge  c > 0;  similarly, 
the  disconnection  of  a module  incurs  an  instantaneous  charge  d > 0. 

While  a module  is  connected,  a rental  charge  of  it  > y^L  per  unit-time 
is  also  incurred. 

Assumption  2.1  is  a restatement  of  the  hypothesis  for  Theorem  1.1. 
Therefore,  using  Theorem  1.1,  the  construction  of  Model  II  will  proceed 
under  the  equivalent  assumption: 

Assumption  2.1'. 

( i)  When  permanent  facilities  satisfy  all  demand,  no  operating  costs 
are  incurred. 

(ii)  When  temporary  facilities  are  necessary  (i.e.,  k < 0)  marginal  costs 
equal  to  the  temporary  facility  charges,  less  an  adjustment 
y^kl  = -y^k  > 0 per  unit-time,  are  incurred. 
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In  (7]  (Example  2.5),  it  was  shown  that  if  the  temporary 


facility  charges  satisfy  Assumption  2.2,  ami  if  the  utilization  of 
temporary  modules  is  required  to  track  the  demand  pattern,  then  Model  I 
of  (7)  was  applicable.  This  situation  can  be  represented  in  a transition 
^ia^ram,  as  shown  in  Figure  2.1.  In  the  diagram,  spares  levels  are  de- 
picted vertically  and  spares  levels  using  the  same,  number  of  temporary 
modules  are  aligned  in  columns.  Each  circled  number  denotes  the  number 
of  temporary  modules  necessary  for  the  corresponding  spares  level  given 
horizontal  to  the  circle  in  the  righthand  margin.  Transitions  occurring 
due  to  departures  are  denoted  by  solid  lines,  while  transitions  occurring 
due  to  arrivals  are  denoted  by  dashed  lines. 

Connection  and  disconnection  (i.e.,  installation  and  removal) 
charges  are  denoted  by  triangles.  Whenever  all  modules  are  fully  used 
(k  * -nL  for  some  integer  n)  and  an  arrival  next  occurs,  two  possi- 
bilities occur.  If  k = -K  (=  -NL  for  some  Integer  N) , then  a permanent 
capacity  expansion  is  undertaken  and  the  N modules  are  disposed  of; 
otherwise  (n  < N) , a new  module  is  connected  at  cost  c,  increasing  the 
total  number  of  modules  in  use  to  n+1.  Similarly,  whenever  a module  is 
being  used  to  serve  a single  customer  (i.e.,  k = -(n-l)L-l  for  some 
n > 1)  and  a departure  next  occurs,  a module  is  disconnected  at  cost  d, 
decreasing  the  number  of  modules  in  use  to  n-1.  In  addition  to  the  in- 
stantaneous connect  and  disconnect  charges,  a cost  rate  of  n per  unit- 
time is  incurred  for  each  module  in  use. 

Given  "free  rein"  on  module  usage,  optimal  disconnections  will 
usually  not  track  the  demand  pattern.  Specifically,  the  optimization 
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dilemma  to  be  resolved  is  whether  or  not  one  should  forestall  a feasible 


disconnection  willfully  incurring  the  penally  of  rental  charges  on  an 
unnecessary  module,  to  protect  against  a subsequent  near-future  reconnection. 
To  illustrate,  consider  Figure  2.2.  Suppose  the  current  spares  level 
is  k = -1;  that  is,  one  unit  of  temporary  capacity  in  a single  module 
is  necessary.  Suppose  that  demand  now  decreases  by  one  (i.e.,  a departure), 
so  that  the  spares  level  passes  to  k = 0;  that  is,  the  temporary  module 
is  no  longer  necessary  to  serve  demand. 

Question:  should  the  module,  now  unnecessary,  be  disconnected 

immediately?  If  the  disconnection  is  made,  then  the  disconnect  charge 
d will  be  incurred  and  there  is  the  risk  of  a near-future  reconnection, 
at  added  expense  c whenever  the  net  demand  (from  now)  next  increases  by 
unity.  On  the  other  hand,  if  the  disconnection  opportunity  is  declined, 
then  the  rental  charges  on  the  unnecessary  module  will  continue  to  accrue. 
Expanding  on  the  above  question,  interest  will  be  focused  on  determining 
precisely  when  module  disconnections  should  be  made. 

With  regard  to  module  connections,  there  is  really  no  question 
involved  with  respect  to  ecoxtomic  optimization.  Assuming  neglible  lead 
times,  it  will  always  be  economically  advantageous  to  defer  connections 
until  necessary  to  serve  demand  (since  c > 0 and  n > )TjL  > 0) . Hence, 
the  connection  path  of  Figure  2.1  is  optimal  and  analysis  will  now  be 
focused  on  the  determination  of  optimal  disconnect  paths. 

2. 2.  Disconnection  Policies 

States  will  be  denoted  by  (n,k),  where  n denotes  the  number 
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of  temporary  modules  currently  connected  and  k denotes  the  current 
spares  level.  Since  all  demand  must  be  served,  (n,k)  is  a feasible 
state  if  and  only  if  nL  > |k|  = -k,  whenever  k < 0.  It  should  be 
emphasized,  however,  that  k is  independent  of  n,  being  only  a function 
of  the  demand  and  permanent  expansion  processes. 

With  regard  to  the  policies  to  be  considered,  note  that  the 
"memoryless"  property  of  the  exponential  distribution  dictates  that 
disconnect  actions  need  only  be  considered  at  demand  epochs.  Furthermore, 
since  the  demand  process  Is  both  time-stationary  and  Markov,  it  suffices 
to  consider  only  time- stationary,  Markov  policies. 

A disconnection  policy  will  in  general  be  denoted  by  u)(’, •)  where 


( 


0, 


u»(n,k)  = J 


if  no  disconnection  occurs  just  prior 
to  entry  into  state  (n,k) 


1,  if  a disconnection  occurs  just  prior 
to  entry  into  state  (n,k). 


Note  that  when  operating  under  a policy  with  u>(n,k)  = 1,  then  state 
(n,k)  will  never,  in  fact,  be  entered,  since  a disconnection  moves  the 
system  to  state  (n-1,  k)  just  prior  to  entry  into  state  (n,k). 


Definition  2.1.  A policy  <«(•,•)  is  feasible  if  and  only  if  cn(0,k)  = 0 
for  all  k and  (n-l)L  > -k  whenever  co(n,k)  = 1,  k < 0.  A feasible 
policy  {!)(•,•)  is  an  s- policy  if  <»>(•,•)  is  determined  by  nonnegative 
scalars  (Sp  Sg,  ...)  such  that 
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w(n,k)  = / 


k < -(n-  1)L  + s , 


k > -(n-l)L  + s , 


n > 1 . 


An  s-policy  is  a constant  s-policy  if  there  exists  a nonnegative  scalar 

s such  that  s = s.  n > 1. 

n ’ — 

Notice  that  Figure  2.1  illustrates  a constant  s-policy  with 

parameter  s = 0;  hence,  this  was  the  form  of  the  disconnection  policy 

considered  in  the  restricted  modular  case  of  Model  I.  Figure  2.3 

illustrates  the  form  of  the  transition  diagram  when  operating  under  a 

constant  s-policy,  assuming  0 < s < L.  Figure  2.4  illustrates  the  form 

of  a constant  s-policy,  assuming  L < s < 2L. 

Physically,  a general  s-policy  can  be  interpreted  as  following 

the  rule:  "whenever  n modules  are  connected,  disconnect  a module  if 

and  only  if  atleast  s units  of  unused  capacity  will  remain  available." 

n 

In  the  next  section,  it  is  shown  that  a constant  s-policy  is  an  optimal 
disconnection  policy. 


2.3.  The  Optimality  of  Constant  s-Policies 

Definition  2.2.  Let  the  random  variable  denote  the  first  time  that 

the  net  demand  reaches  i: 

= min(  t > 0 : fl(  t)  - D(0)  = i},  i = 0,  _+l,  +2f  ...  . 

In  the  absence  of  intervening  permanent  expansion,  the  excess  demand 
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(-k)  tracks  demand.  Hence,  if  there  ts  no  intervening  permanent  expan- 
sion, T.  is  also  the  first  time  that  the  spares  level  decreases  by  an 
amount  i,  i > 0.  Similarly,  if  there  is  no  intervening  permanent  expan- 


sion, T is  the  first  time  that  the  spares  level  increases  by  an 
amount  i i > 0. 

Definition  2.1.  Define  h(*,*)  hy 

-rT 

h( u s)  = E[e  U|l  < T ] P(T  < T } 

* ' u s u s 

Define  H(*,*)  by 

( 0,  u = s 

+■  ( c + -p)  h(  u+1,  u-s)  - (d  - ■— ) h( u-s,u+l) , s > u > 0 

To  understand  the  role  of  H( • , • ) with  regard  to  disconnections, 
suppose  that  the  system  is  in  state  (1  u- 1) , u > 0.  and  that  a departure 
now  occurs.  It  is  desired  to  economically  compare  two  alternative  courses 
of  action: 

Action  1 ; immediately  disconnect  the  unnecessary  module;  or 

Action  2:  do  not  disconnect  the  unnecessary  module  until  a departure 

occurs  when  the  system  is  ii.  state  (1,  s-1),  s > u. 

Figure  2.5  depicts  these  two  alternative  actions.  As  illustrated  there, 
the  system  will  be  "restored"  to  the  same  state  whenever  either  of  the 
states  (0,s)  or  (1,-1)  are  again  first  reached  --  these  two  states 
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will  be  called  the  reentry  states.  The  corresponding  incremental  reentry 


time  is  obviously  distributed  as  min(T  , . T ).  To  make  the  economic 

u+1*  u-s 

comparison,  the  expected  savings  which  accrue  from  following  Action  2 

instead  of  Action  1 will  be  computed  over  this  reentry  interval  of  time. 

The  calculation  is  done  by  conditioning  on  the  reentry  state. 

If  T^  < T^  then  state  fO,s)  is  the  reentry  state.  Following 

Action  1 an  immediate  disconnection  charge  d is  incurred;  following 

Action  2,  the  rental  charge  n per  unit- time  is  incurred  over  the 

interval  [0.  T J and  a deferred  disconnection  charge  d is  incurred 
* u-  s 

at  T^  . Hence,  the  savings  of  Action  2 over  Action  1 are  given  by 
T 

u-s 

S = d - (/  n e'rt  dt  + d e'rTu-8) 

0 

= (d  - - (d  - -)  e’rTu_s  if  T < T . (2.1) 

v'  ' r'  * u-s  u+1  ' 


If  T , < T , then  state  (1.  -1)  is  the  reentry  state. 

Following  Action  1 an  immediate  disconnection  charge  d is  incurred 

and  subsequently,  a reconnection  charge  c is  incurred  a the  reentry 

time  Ty  Following  Action  2,  only  the  rental  charge  it  per  unit-time 

is  incurred  over  the  time  interval  [0.  T Thus,  the  savings  of 

’u+1  ’ 

Action  2 over  Action  1 are  given  by 


-rT,  , u+1 

S a (d  + ce  U+  ) - / n e'r  dt 

0 


-rT, 


- (d  -i)  + (c  +i)  e'“u*1  , 


if  T . < T 

u+1  u-s 


(2.2) 
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Using  (2.1)  and  (2.2)  the  expected  savings  of  Action  2 over 


Action  1 are  given  by 

ES  = E[S|T  < T ] P{T  < T ,}  + E[S)T  . < T ] P(T  , < T ) 
u-  s u+ 1 u-s  u+ 1 u+  1 u-s  1 u+ 1 u-  s ; 

= EKd  - 7)  - (d  - 7)  e Tu"s|t  < T J P(T  < T ,) 
r r u-s  u+1  u-s  u+1 

+ E[(d  - ^)  + (c  + 7)  e"rTu+l|T  < T ] P(T  < T } 
r r u+1  u-s  u+1  u-s 

It  . < T ] P(T  , < T ) 
u+1  u-s  u+1  u-s' 


•(<>-?)♦(«♦*)  E[e-W' 


-rT 


- (d  - 7)  E[e  U"S|T  < T ,]  P{T  < T .) 
r u-s  u+1  u-s  u+1 

= (d  - “)  + (c  +7)  h(u-*-l,  u-s)  - (d  - •£)  h(u-s,  u+1) 


= H(u,s) 


(2.3) 


Tlius^  H(u,  s)  represents  the  expected  savings  accruing  from  a single 
disconnection  deferral  up  to  spares  level  s when  an  opportunity  to 
disconnect  into  state  (0,u)  is  at  hand.  In  particular,  the  expected 
savings  accruing  from  the  disconnection  deferral  previously  illustrated 
in  Figure  2.3  is  given  by  H(0,  s). 


Lemma  2. 1. 

H(u,s)  = H(u,  v)  + H(v,s)  h(u-v,  u+1),  0 < u < v < s . 

Proof.  The  reentry  states  are  (0,s)  and  (1,-1),  with  respective 

reentry  states  T and  T 

u-s  u+1 
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Case  1:  T < T In  this  case,  spares  level  v is  reached  prior 

u-v  u+1  ’ 

to  reentry  since  v < s implies  T < T Since  0 < u < v < s, 

r u-v  u-s  — • 


T ? T + T' 
u-s  u-v  v-s 


and 


u+1 


St 


u-v 


T'  . 
v+1 


> 


where  (T|)  is  a sequence  of  random  variables  independent  of  (T^)  and 
with  distributions  identical  to  those  of  fT.l. 

i 


Case  la:  T'  < T'  T < T .).  In  this  subcase,  reentry  in 

v-s  v+1  v u-s  u+1'  ’ 

state  (O.s)  occurs  at  reentry  time  T + T'  Proceeding  as  before, 

we  have 


S = 


■ rTv- s ~ r^u- v 


> 


if 


T 

u-v 


< T 


u+1 


and 


T' 

v-s 


< T' 
v+ 


l * 


(2.4) 


Case  lb:  T*  < T'  . (=>  T , < T ).  In  this  subcase,  reentry  into 

v+1  v-s  ' u+1  u-s'  ’ 

state  (1,  -1)  occurs  at  reentry  time  Tu  v + T'  Proceeding  as  before, 
we  have 


c , . Jtx  / ju  "rTv+l  ”rTu-v 
S=(d-“)+(c+-)e  e 


if  T < I , and  T'  . < T1 

U-V  u+1  v+1  V-  3 


(2.5) 
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jL 


Case  2:  T 


. < T . 
u+1  u-v 


at  reentry  time 
Hence,  as  before, 


u+1' 


In  this  case,  reentry  into  state  ( 1,  -I)  occurs 
without  an  intermediate  visit  to  spares  level  v. 


if  T 


u+1 


u-v 


(2.6) 


For  notational  convenience,  denote  the  pertinent  events  by 

£.  = {T  < T = (T'  < T'  . ),  = {T*  . < T'  ),  and 

1 u-v  u+1  9 la  v-s  Vfl  9 lb  v v+1  v-s  9 

e2  = (Tu+l  < Tu-v}‘  Note  that  P(£la  U &!1  £1}  = P<£la]'  P{£1  U £lb|£l] 
= F(£lb)  and  P{£la)  + P^£lb^  = Also  denote  a = d - jt/r  and 

b = c + n/r.  Rewritting  (2.4)  - (2.6)  in  new  notation  gives 


- rT  1 - rT 

tJ-v-s  LAu-v 


/ a - a e 


e / K *rTV+l  *rTU-V 

S = ( a + b e e 


a + b e 


rTu+l 


if  ex  u ela 


if  £]L  u elb 


if 


h. 


Thus, 
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H(u,s)  = ES 


E[a  - a e lTv"S  e’rTu_vlei  U Elfl]  PUj  U £la) 


+ E [a  -t-  b e tTv+1  e rTu_V|g1  U &lbJ  P{  U e]b) 


+ E[a  + b e U+^lgg]  Ple^j 


r E[a  - a e'rT"-S  e’rTu'V|gla]  P^j 


= E 


I-  E[a  + b e‘rT;+1  e‘rTu'V(glb]  Pf^] 


Pf^} 


+ E[a  -i-  b 


frTuale2]  Pt^} 


= E 


E[a  - a e 


-rT 1 


v-s  | 


eul  nela) 


• rT, 


E[a  + b e“rTv+1|&lb]  P{eib) 


+ (a  - a e'rTll-v)  (P{eia}  + PUlb}) 


Pl^J 


+ E[a  + b e'rTu+1Iep]  P{ &2 ) 


e[h(v,s)  e'rI’u'v!e1]  p(e1) 


E [a 


a e 


'rTu-v|e1]  P{e1)  + E[a  + b e"1 Tu+l | Ep]  P{g  2) 


= H(Vj3)  h(u-v?  u+1)  +H(u>v) 
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Definition  2.4.  Let 


s*  denote  the  greatest  supremum  of  H(0,s)  : 

s*  = sup{s'  : 11(0,8')  = sup  H(0,s)J. 

s > 0 

Lemma  2.2.  s*  < + oo  d < it/ r. 

Proof.  The  proof  of  necessity  (=>)  is  by  contradiction.  Suppose  d >n/r. 
Then,  for  s > 1, 

H(s-1,  3)  = (d  - £)  + (c  + i)  h(s,-l)  - (d  - i)  h(-l,s) 

= (d  - f)  (1  - h( - 1,  s))  + (c  +i)  h(s,  -1)  >0  , 

since,  by  Definition  2.3,  0 < h(u,v)  < 1 for  u < 0 < v and  v < 0 < u. 
Thus,  using  Lemma  2.1 

H( 0, s)  = H( 0,  s-1)  + H( s- 1,  s)  h( 1-s,  1)  > H(0,  s- 1) , s > 1, 

since  0 < h(l-s,  1)  < 1 for  s > 1.  Hence,  d > =»  s*  a contra- 

diction. Thus,  the  original  supposition  was  false  and  we  must  conclude 

To  prove  sufficiency  (<^) , suppose  d < Then 

H(s-1,  .)  - (d  - f)  (1  - h(-l,s))  + (c  + £)  h(s,  -1) 

= (d  -f)  (1  - asPs)  + (c  +f)  bs(l-Ps)  , (2.7) 
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- rT  i 

where  ag  = E[e  " |l  i<  Tg],  bg  = E f e~ rTs | Ts  < T L , and  Pg  = P(T  l <Ts). 

Since  A.p.q  > 0.  it  follows  that  a decreases.  P increases  and  b 
* * * s ' s s 

decreases.  Furthermore.  P -*  1 as  s ->  +<».  Since  0 < b <1 

* s s’ 

b (1-P  ) — > 0 as  s ->+<».  As  (c  + — ) b ( 1— P ) ->0,  J s'  such  that 

s s'  ' r'  s s'  ’ 

0 < (c  + ~)  bg(  1-Pg)  < (“  - d)  ( 1-^),  for  aH  s > s *.  Hence,  from 

(2.7),  H(s-1,  s)  < 0,  for  all  s > s'.  Thus,  using  Lemma  2.1,  we  conclude 

that  H(0,s)  < H(0,  s- 1) , for  all  s > s'  and  it  follows  that  s*  < s ' < +00. 

□ 

The  economic  interpretation  of  Lemma  2.2  is  simple.  The  lemma 
states  that  disconnection  will  be  advantageous  if  and  only  if  the  dis- 
connection charge  is  less  than  the  cost  of  renting  a module  over  an 
infinite  horizon.  Henceforth,  it  will  always  be  assumed  that  d < — . 

The  next  lemma  will  prove  to  be  quite  important  from  a computational 
standpoint  (see  Section  2.8). 

Lemma  2,3.  s*  is  independent  of  X,  K,  N,  g(*). 

Proof.  Since  T ^ and  Tg  (s  > 0)  are  dependent  only  on  the  demand 
process,  H(0, s)  is  completely  determined  by  the  demand  process  and  the 
parameters  c,  d,  jt  and  r.  □ 

Leircna  2.4.  H(s,s*)  >0,  s < s*  and  H(s*,s)  < 0,  s > s*. 
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Proof.  In  view  of  Lemma  2.1  and  Definition  2.3(  anything  to  the  contrary 


would  contradict  Definition  2.4. 


□ 


The  next  two  lemmas  establish  bounds  on  the  behavior  of  h(  • • ) 
with  regard  to  each  of  its  arguments. 

Lemma  2.5. 


-rT. 


h(u,s.l)  . h(u,s)  ♦ Efe  U|TS  < T,  < Tg+1 ] P(Tg  < Tq  < Xs+1)  , 


u < 0,  s > 0. 


Proof. 


h(u,s+l)  . E[c'rT‘‘llu  < TstlJ  mu<Ta+1) 


-rT 

= E[e  UlT  < T . T < T ,]  P{T  < T . T < T .} 

1 u s’  u s+1  1 u s>  u s+1 

— rT 

+ E[e  U|T  <T  <T  J P(T  <T  <T  ,} 
s u s+1  s u s+1 

- rT 

= E [e  U[T  < T ] P{T  < T } 

' u s u sJ 

- rT 

.Ele  “lTs  < Iu  < Tstll  P(I,  < T„  - t>+1) 

* h<“,s>  + E[c'rT“|lt  < Iu  < T t]  P(TS  < Tu  < Ts+1) 


The  next- to- last  equality  follows  since  {T^  < Tg}  =i>  (Tu  < Tg^}.  □ 


3J 


Lemma  2.6. 


-rT, 


h(u,s)  > h(u+l,s)  + E[e  U|T  <T  <T  , ] P{T  <T  <T 

' ’ ' u s u+ 1 u s u+ 1 ’ 

u > 0;  s < 0 . 


Proof. 


-rT 


h( u s)  = E[e"  UjT  < T ] P{T  < T ) 
’ ' u s u s 


■rT, 


= E[e  U|T  < T , T . < T ] P(T  < T , T , < T } 
u s'  u+1  s u s’  u+1  s 


-rT, 


E [e  U|T  < T < T J P{T  < T < T ,} 
u s U+-1  u s u+lJ 


-rT,, 


= E[e  U|t  . < T ] P{T  , < T ) 

U+1  s u+1  s 

- rT 

+ E [e  |t  < T < T , ] P{T  < T < T 

u s u+1  u s u+lJ 

> Efe  rTu+1|T  < T ] P{T  . < T ) 
u+1  s u+1  s 

-rT, 


uIt 

< T 

< x , 

u 

s 

u+1 

-rT. 

U 

< 

T < T 

u s u+1  u s u+lJ 


The  third  equality  follows 

from 

(T  , 
u+1 

< T } =>  (T  < T } for 
s u s 

u > 0. 

_ 

rTu  - rT, 

The  inequality  follows  since: 

(i)  for 

u > 0.  T . > T =>  e 
— ’ u+ 1 u 

->  e 

(with  probability  1);  and 

(ii) 

P(T  . 
u+1 

< T } > 0 for  u > 0. 
s — ’ 

s < 0.  □ 

The  economic  interpretation  of  the  following  theorem  is  that  the 
marginal  savings  of  disconnection  deferrals  are  decreasing. 
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Theorem  2.7.  H(s,  s+1)  is  decreasing,  s > 0. 


Proof. 

H(s,  s+1)  = (d  - + (c  + 7)  h( s+1,  -1)  - (d  - ^)  h(-l,  s+1)  , 

H(®-1,  s)  = (d  - f)  + (c  + f)  h(s,  -1)  - (d  - i)  h(-l,  s)  . 

Thus, 

H(s,  s+1)  - H(  s- 1,  s)  = (c  +7)  (h(  s+1,  -1)  - h(s,  -1)) 

- (d  - 7)  (h(-l,  s+1)  - h(-l,  s))  . 

To  prove  the  theorem,  it  suffices  to  show  that  the  right-hand  side  of 
the  above  expression  is  negative.  Using  Lemmas  2.6  and  2.5  gives 

-rT 

h(  s,-l)  - h(s+l  -1)  > Efe  SIt  < T < T .]  P{T  < T , < T ) 

s -I  s+1  s -1  s+1 

> E(*'rT-1iT3  < I,  < Ts+1]  P{Is  < < T>+1 

= h(-l,  s+1)  - h(-l,  s)  . 

Also,  (c  +7)  > -(d  - 7)  = (7  - d)  > 0.  Combining  this  with  the  above 
inequality  yields 


(c  + 7)  (h(s,  -1)  - h(s+l,  -1))  >-  (d  - i)  (h(-l,  s+1)  - h(-l,  s)) 


which  implies 

(c  + (h(s+l,  -1)  - h(s-l))  - (d  - -p)  (h(-i,  s+l)  - h(-l,  s))  <0  . 

□ 

Definition  2.^.  H(0,*)  will  be  said  to  be  strictly  unimodal  about 

[s',  s"]  if:  (i)  11(0, s)  < H(0,  s+l),  s < s ( ii)  H(0,s)  <H(0,  s+l), 

s'  < s < s";  and  (iii)  H(0,s+1)  < H(0,s),  s > s". 

The  following  theorem  will  be  important  later  with  regard  to 
computing  s* . 

Theorem  2.8.  H(0,*)  is  strictly  unimodal  about  [s^-1,  s*J. 

Proof.  From  Lemma  2.4,  H(s*-1,  s*)  > 0 > H( s*,  s*+l) . Therefore,  by 
Theorem  2.7,  H(s-1,  s)  > O for  s < s*  and  H(s,  s+l)  < 0 for  s > s*. 

Hence,  using  Lemma  2.1  and  the  fact  that  h(-s,l)  > 0 for  all  s > 0, 

/ | 
> H(0,s),  s < s*-l 

H(0,  s+l)  =-  H(  0,  s)  + H(s,  s+l)  h(  -s,  1)  />H(0,  s*- 1) , s = s*-l 

< H(0,s),  s > s*  □ 

Thus,  the  maximizing  set  of  11(0, •)  contains  either  the  single 
component  s*  or  the  two  components  s*-l  and  s*.  The  optimal  of 
the  constant  s-policy  with  parameter  s*  will  now  be  established. 
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Lemma  2.Q.  If  u)  is  a feasible  disconnection  policy  and  if,  for  any 
n > 1,  inf(i  : u>(n,  -(n-l)L  + i)  = 1)  > s*,  then  ou  is  not  optimal. 

Proof.  Suppose  u>  and  n are  such  that 


s = inf{i  : to(n,  -(n-l)L  + i)  = 1)  > s*  , 


To  demonstrate  the  non-optimality  of  to.  we  construct  an  alternative 

policy  cn'  (possibly  nonstationary)  realizing  expected  costs  strictly 

less  than  those  of  co.  Specifically,  define  a>'  as  follows: 

/ 

0,  k < -(n-l)L  + s* 

1,  k = -(n-l)L  + s* 


(a)  o>'(n,k)  = 


( b)  Whenever  in'  disconnects  from  (n,  -(n-l)L  +■  s*-l)  to 

(n-1,  -(n-l)L  +•  s*),  a'  makes  no  further  disconnections  until 
(just  prior  to)  first  reaching  either  (reentry)  state 
(n,  -(n-l)L-l)  or  ( reentry)  state  (n-1,  -(n-l)L+s). 

(c)  o'  is  otherwise  identical  to  co. 

If  s = then  (reentry)  state  (n-1,  -(n-l)L  + s)  is  interpreted  as 
an  unreachable  state  in  (b).  If  s < +o°  and  u>(n-l,  -(n-l)L+s)  = 1, 
then  (b)  and  (c)  are  interpreted  to  mean  that  to '(n-1,  -(n-l)L+s)  = 1. 

Because  of  (b),  oo'  need  not  be  stationary.  However,  we  can  still 
represent  the  two  alternative  strategies  in  a transition  diagram  by 
defining  "dummy"  states  {(n-l,k)',  -(n-l)L  < k < -(n-l)L  + sj,  asshown 
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in  Figure  2.6.  Given  this  enlargement  of  the  state  space,  o>  and  ca' 
be  represented  by 


can 


cn(  n-  1,  ( n-  1)L  + s) 
o)'(n-l,  k)  ' 

and 

o>'(  i,k) 


= 1 , 

0,  -(n-l)L  < k < -(n-l)L  + s 

" < 

1,  k - -(n-l)L  + s 
i / n 

i = n,  k < -(n-l)L  + s* 

i = n,  k = -(n-l)L  + s* 

i = n,  k > -(n-i)L  + s*  , 


where  it  is  understood  that  a>'  disconnects  from  ( n,  -(n-l)L  + s*  - 1) 
to  (n-1,  -(n-l)L  + s*) 

To  compare  a>  and  a>'  we  suppose  that  the  system  has  been  in 

state  (n,  -(n-l)L  + s*  - 1)  and  that  a demand  decrease  now  occurs.  As 

illustrated  in  Figure  2.6,  policy  id'  immediately  disconnects,  moving 

to  state  (n-1,  -(n-l)L  + s*)  ' . Under  policy  u>,  the  disconnection  is 

deferred  until  s-s*  additional  demand  decreases  occur.  At  time 

T - min{T_^  , . T „ }.  the  system  will  be  "restored"  to  identical  condi- 

s*-s  * 

tion  regardless  of  which  strategy  is  chosen.  Hence,  any  expected  savings 
of  oj  over  o>'  must  accrue  over  the  time  intervals  of  the  form  [0,T] 
(since  u>'  is  identical  to  o>  elsewhere). 
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If  T then  policy  to  incurs  an  additional  rental  charge 

of  it  per  unit-time  over  the  interval  [0,T],  while  policy  to'  incurs 
an  immediate  disconnexion  charge  d,  followed  by  a reconnection  charge 
c at  T . (Notice  that  the  adjustments  — y _ | k I do  not  enter  in  the 
dif ferential,  since  the  spares  level  is  not  altered  by  the  choice  of 
disconnection  policy;  both  policies  will  incur  equal  adjustments.)  Hence, 


(d 


■rT 


+ ce 


S*4-l 


S*+l 


) - 


ne 


-rt 


dt 


(d  - f)  + (c  + 7) 


- tTs*+  1 


if  T , < T 

s*+l  s*-s 


(2.3) 


If  T = T (s  < e<»),  then  policy  cn  incurs  an  additional  rental 
charge  of  it  per  unit- time  over  the  interval  [0,T]  followed  by  a dis- 
connection charge  d at  T . Policy  to'  incurs  only  an  immediate 
disconnection  charge  d.  Hence, 

fTs*-s  _rt  -sTs*.s 

S s d - (/  ice  dt+de  ) 

0 


- (d  - f)  - (d  - -J) 


-rT 


s*-s 


if  T 


s*-s 


< T 


s*+ 1 


(2.9) 


(2.8)  is  (2.2)  evaluated  at  u = s*  and  (2.9)  is  (2.1) 
u s*.  Hence,  from  (2.3),  the  expected  savings  of  to 
given  by  ES  = H(s*,s).  But,  by  Lenina  2.4,  H(s*,s)  < 0 
Furthermore, 


evaluated  at 
over  to'  are 
for  s*  < s < 
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H(s*f»)  - (d  - i)  + (c  + i)  E[ 


-rTS*fl  ; 


since  (2.9)  is  positive  for  all  s > s* ; hence,  H(s*,  00)  < 0.  Thus,  a) 
realizes  (strictly)  negative  expected  savings  over  a>';  that  is,  to' 
realizes  expected  costs  strictly  less  than  those  of  cl>.  Hence  cu  is 
not  optimal.  □ 


Lemma  2. 10.  There  exists  an  optimal  disconnection  policy  a#  with 


co*(  l,k)  = 


k < s* 


k = s* 


Proof.  Referring  to  Figure  2.2,  suppose  that  the  system  has  been  in 
state  (1,  -1)  and  that  a departure  now  occurs  (i.e.,  k ->  0 from  -1). 
Compare  the  two  alternative  strategies  of  either  an  immediate  disconnection 
or  a disconnection  deferral  until  s > 0 units  of  spare  capacity  will 
remain.  The  savings  accruing  from  the  deferral  strategy  is  H(0,s) 
s > 0.  We  can  do  no  better  than  maximize  this  savings  with  H(0,s*), 
by  assigning  co*(l,*)  as  above.  □ 


Theorem  2.11.  There  exists  an  optimal  disconnection  policy  u#  with 


oj*(n,k)  = 


k < -( n-  1)L  + s* 
k = -(n-l)L  + s* 
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Proof  (Induction).  By  Lemma  2.10,  the  theorem  is  true  for  n = l and 


we  now  assume  that,  in  general,  it  is  true  over  1,  2,  ...,  n-1,  for 
some  n- 1 > 1 . 

Let  <ij  be  an  arbitrary  policy  agreeing  with  (2.10)  over 
1,  2,  ...,  n-1.  Let 

s = inf(i  : co(n,  -(n-l)L  + i)  = 1)  . 

By  Lemma  2.9,  it  suffices  to  assume  that  s < s*.  Let  a,'  be  a policy 

agreeing  with  u>  everywhere  except  for  io'(n,  -(n-l)L  + i)  = 0,  0 < i < s* 

and  a>'(n,  -(n-l)L  + s*)  1.  We  shall  show  that  the  expected  savings 

of  cjd  1 over  an  are  nonnegative.  If  s = s*,  then  a>'  s :o  and  there 

is  nothing  to  prove;  hence  assume  s < s*. 

The  two  alternatives  are  illustrated  in  Figure  2.7.  Suppose,  the 

system  is  in  state  (n,  -(n-l)L  + s - 1)  and  that  a demand  decrease  now 

occurs.  Let  T = minfT  ..  T Under  policy  m an  immediate  dis- 

s+l'  s-s*  ’ 

connection  to  state  (n-1,  -(n-l)L  + s)  occurs.  Furthermore,  under 

policy  o>,  the  system  actual  enters  state  (n-1,  -(n-l)L  s)  and  no 

further  disconnections  are  made  over  the  time  interval  [0,TJ  since,  by 

the  induction  hypothesis,  j>(n-l,  (n-l)L  + i)  --  u>(n-l,  (n-2)L  + (i-L))  = 0 

for  i < L+s*  (and  s*  < L+s*) . Under  policy  a'1,  the  disconnection  is 

deferred  until  s*  units  of  spare  capacity  will  remain.  At  time  T, 

the  system  is  restored  at  one  of  the  two  reentry  states  (n,  -(n-l)L  - l) 

or  (n-1,  (n-l)L  + s*)  (regardless  of  which  policy  is  chosen)  depending 

upon  whether  T=T  , or  T^T  respectively. 

r s+1  s-s*> 
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If  T = Tg  then  a>  incurs  an  immediate  disconnection  charge 

d and  a subsequent  reconnection  charge  c at  time  T u)'  incurs 

s+1 

an  addition  rental  charge  of  « per  unit-time  over  the  interval  [0,TJ. 
(Notice  that  the  adjustments  -r^lkl  do  not  enter  in  the  differential, 
since  the  spares  level  is  not  altered  by  choice  of  disconnection  policy; 
both  policies  will  incur  equal  adjustments.)  Hence, 

T i 

-rT  i s+1 

„ ,,  As+L  , -rt 

S = ( d + ce  ) - / ne  dt 

0 


(d  -f) 


/ itv  -rTs+l 
+ ( c + -)  e 


if  T . < T „ 
s+1  s-s* 


(2.U) 


If  T = 
charge  d.  cjd' 
the  interval 


T . Hence 
s-s*  > 


T then  to  incurs  only  the  immediate  disconnection 

s-s*> 

incurs  the  rental  differential  n per  unit- time  over 
[0,  Tg  g^j  followed  by  a disconnection  charge  d at 


S = d 


(/ 

0 


s-  s* 


ne 


■ rt 


dt 


de""18-^) 


= (d  - f)  - (d  - “)  e 


■r^s-s* 


if  T „ < T . 
s-s*  s+1 


(2.12) 


Comparing  (2.12)  with  (2.1)  and  (2.1l)  with  (2.2),  it  follows  from  (2.3) 
that  ES  = H(  s,  s*).  But,  by  Lemma  2.4,  H(  s,  s*)  > 0,  since  s < s*. 
Hence,  nonnegative  expected  savings  accrue  from  using  co'  over  00. 

Since  <0  was  arbitrary,  it  follows  that  the  theorem  holds  over 
1,  2,  , n.  The  theorem  then  follows  by  induction  on  n.  □ 

Thus,  the  optimality  of  constant  s-policies  is  established. 
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2.4.  The  Transition  Equations  for  Constant  s- Policies,  s > 1 

Recall  that  K denotes  an  upper  limit  on  the  temporary  facility 
capacity  that  can  be  used.  Since  the  temporary  facilities  are  modular 
of  size  L,  only  those  values  of  K being  nonnegative  integer  multiples  of  L 
are  of  interest  here,  Hence,  K will  be  denoted  as  K = NL,  where  N 
represents  the  limiting  number  of  temporary  modules  allowed. 

cj^X,  N,  s)  will  denote  the  expected  total  discounted  cost  when 
the  initial  spares  level  is  k,  the  initial  number  of  temporary  modules 
connected  is  n,  the  limit  on  temporary  module  usage  is  N,  the  permanent 
capacity  expansion  size  is  X+l,  and  when  the  module  disconnection  pattern 
adhers  to  the  constant  s-policy  having  parameter  s;  X > NL,  N > 0, 
s > 1,  k > -NL,  0 < n < N.  The  constant  s-policy  having  parameter  s = 0 
is  not  considered  here,  since  that  case  has  previously  been  treated  by 
Model  I of  [7]. 

Similarly,  F^  will  denote  the  expected  incremental  discounted 
costs  (until  the  next  arrival  or  departure)  while  the  spares  level  has 
value  k,  when  there  are  n modules  connected.  Notice  that,  by  Assumption 
2.2,  these  costs  are  independent  of  N.  Also,  by  Assumption  2.1’, 

F°  = 0 for  k > 0 ard  F^  includes  an  (expected)  adjustment  term 

K K 

oo  - - - 

/ r,kr'  (1  - e”rt)A  e"  1 dt  = (A+r)"  r,k  , k < 0 . 

0 

Recall  that  p is  the  probability  of  arrival  and  q = 1-p  is 
the  probability  of  departure.  Denote: 
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O.  = Aq(  A+r) 


and  p = Ap(A+r) 


(2.13) 


-1 


Also,  denote: 

II  = (A+r)'1  it  and  r = (A+r)"1  . (2.14) 

When  a constant  s-policy  is  in  use,  there  are  four  different  types 
of  states:  (a)  potential  module  connection  states,  ( b)  potential  module 

disconnection  states,  ( c)  states  where  neither  connection  nor  disconnection 
possibilities  exist,  and  (d)  a potential  permanent  expansion  state.  The 
incremental  expected  costs  for  each  of  these  state  types  are  derived  below. 

(a)  Potential  connection  states  are  of  the  form  (n,k)  where 

k = -nL,  0 < n < N-l.  At  these  states,  an  arrival  triggers  the  connection 
of  an  additional  module  at  cost  c.  Conditioning  on  the  time  and  the 
type  of  the  next  event  gives 

00 

- / (pc  e n + (nit  + ky,)  r‘  L(  1 - e*rt))X  e"^1  dt 
0 1 

= ^c  + nil  + kl  , k = -nL,  0 < n < N-l 

( b)  Potential  disconnection  states  are  of  the  form  (n,k)  where 
k = -(n-l)L  + s - 1,  1 < n < N.  At  these  states,  a departure  triggers 
the  disconnection  of  a module  at  cost  d.  Notice  that  these  states  are 
distinct  from  those  of  (a),  since  s > l.  Conditioning  on  the  time  and 
the  type  of  the  next  event  gives 


Fk  = ^ e + (nn  + min(0,k)  y^)  r k(  1 - e rL))A  e~ ,'t  dt 

= Od  + nil  + min(0,k)r  , k = -( n- l)L  + s - 1,  1 < n < 

(c)  Intermediate  states  where  there  is  no  potential  for  either 
connection  or  disconnection  are  of  the  form  (n,k),  where  -nL+1  < k 
< -( n- l)L+s-2,  1 < n < N.  These  states  are  distinct  from  those  of  (a) 
and  (b).  However,  in  the  special  case  of  s = L = 1,  these  states  will 
be  nonexistent.  Conditioning  on  the  time  of  the  next  event  gives 


= nil  + min(0,k)l  , -nL+1  < k < -(n-l)L+s-2,  1 < n < N . 

(d)  The  potential  expansion  state  is  (N,  -NL)  and  an  arrival 
in  this  case  triggers  an  expansion.  Conditioning  on  the  time  of  the 
next  event  gives 


/ (N.t  - NL  y ) r_1(l  - e"rt)A  e dt 
0 1 

* Nn  - NLr 

When  expansion  occurs,  there  may  be  extraordinary  costs  ( seperate  from 
the  expansion  cost)  for  disposing  of  the  temporary  facilities.  These 
costs  will  be  denoted  0^(N).  For  computational  reasons,  it  will  prove 


convenient  to  write 


f-nl  = Pc  + Nl[  ' NLr 


It  will  be  demonstrated  sho-lly  that  this  causes  no  inaccuracies  provided 
that  a term  "-c"  is  included  in  the  disposal  cost  function  0 (•).  For 
Example,  0g(N)  !«d-c  would  account  for  the  disconnection  of  the  N 

modules  when  permanent  expansion  occurs,  as  well  as  rectifying  the  non- 
existent connection  charge  included  in  the  above  expression  of 

-NL 

It  will  also  be  assumed  tlut  is  nondecreasing  and  that  c + 0 (.) 

is  nonnegative  (but  otherwise  arbitrarily  given). 

Summarizing,  the  expected  incremental  costs  are  given  by: 


/°< 


>!- 


k > 0,  n = 0 


k = -nL,  n > 0 


fJc  + nil  v kr  , 

+ nil  + min(0,k)r  , k = -( n- l)L+-s- 1,  n > 1 


' nil  + min(0,k)l  , 


-nL+1  < k < -(n-  1) L+s-2,  n > 1 


By  conditioning  on  the  time  and  type  of  the  next  event,  the  cost 
transition  equations  have  the  general  form 


00 

c£(X,N,s)  = f£  + / ^^^(1,11,5)  + pC^X^s))  e'rt  dt 


= Fk  + ^l^A*)  + PCj"i(X,N,i)  , 
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where  n'  t (n-i  n)  and  n"  « (n,nfl)  depending  on  whether  or  not 
either  a connection  or  a disconnection  state  is  involved.  Given  the 
expected  incremental  costs  derived  in  (a)  - (c),  the  cost  transition 
equations  for  a constant  s-policy  are: 

C°(X,N,s)  = QC°+1(X,N,s)  + PC^X^s),  k>0  (2.15) 

<(X,N,S)  = F"  + «cJ+i(x,N,s)  + f3C^JfX,N,s), 

k = -nL,  0 < n < N-l  (2. 16) 

c£(X,N,s)  = Fr  + QC^^X»N»fl)  + PCk.l(X’N«S^ 

k = -( n- l)L+s- 1,  1 < n < N (2.17) 

C^(x,N,s)  = f"  + qc"+1(x,n,s)  + 

-nL+1  < k < -(n-l)L+s-l,  1 < n < N . (2.18) 

Notice  that  no  upper  bound  is  imposed  upon  k in  (2.15).  This  corresponds 
to  the  "untruncated  demand  assumption"  discussed  in  Chapter  1.  Therefore, 
the  comments  in  [7]  (Section  1.7)  regarding  the  viability  of  this  assumption 
are  applicable  here. 

When  k = -NL  and  an  arrival  next  occurs,  a permanent  expansion 
of  size  X+l  (X  > NL)  is  undertaken  at  cost  g(X).  Whenever  this 
happens,  the  spares  level  increases  to  the  value  X-NL  > 0 and  all 
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temporary  modules  are  disposed  of  at  cost  p^(N).  Therefore,  conditioning 
on  the  next  event  time  and  type,  the  boundary  transition  equation  becomes 

OM'-o 

CD 

^ + f Ae'M^NNL,l(X>N'3)  + p(Pe(N)  f 8(X)  + CX-NL(X^N»S)))e'rt  dt 

pNNL  + QC-NL+1(X»N'S)  + P^e(N)  + 8(X)  + CX-NL(X’N'S)]  ' (2'193 

Notice  that  the  inclusion  of  the  "-c"  term  in  (N)  will  indeed  cancel 

N 

the  nonexistent  "(3c"  charge  term  of  F . 

The  solutions  to  equations  (2.15)  - (2.I9)  are  treated  in  the 
next  four  sections.  A final  note:  when  N - 0,  the  system  of  equations 

consists  of  only  (2.15)  and  (2.19). 
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2.5.  Solutions  £j>r_the_Trans  i t ion  K<Luatlj.ms,  k '■>  0,  n = 0 

The  special  nature  of  the  linear  equations  (2.15)  are  well-known 
[1],  [3].  As  the  solutions  to  these  equations  are  derived  in  (7),  the 
corresponding  results  are  simply  listed  here. 

Lemma  2.12.  For  a ani  [1  given  by  (2.13): 

( i)  a > 0,  (3  > 0, 

(ii)  < l, 

( iii)  Op  < 1/4. 

Proof.  See  Lemma  2.1  of  [7). 

Theorem  2.13. 

C°(X,N,s)  = Zk  C°(X,N,s)  , k>0, 

where 


2 

Furthermore,  Z is  real  and  |3  < Z < Q+p.  Also,  QZ  - Z + (3  _ 0. 
Proof.  See  Theorem  2.2  of  [7].  O 


(2.20) 


(2.20a) 
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Corollary  2. I ..  For  X > NL,  X ' > N 'L,  s > 1,  s ' > 1 and  k > 0, 


co  x>N»s)  <co(x,>N,»s’)^c^x»N^)  <<£(X',N',s')  . 


• 7,10  result  is  a direct  consequence  of  the  theorem  since  Z > p > 0.  □ 

As  shown  in  [7],  Z1  Is  the  Laplace  transform  of  , evaluated 
at  the  discount  rate  r,  for  i > 1. 

i -rTi 

Lemma  2,1^.  Z = E[e  ],  i > 1. 


Proof.  See  Lemmas  2.4  and  2.5  of  [7]. 


2.6.  The  Form  of  the  Functionals  {c”(*,N, s)} 

When  temporary  modules  are  in  use,  the  expected  costs  behave 
according  to  (2.16)  - ( 2 . I9) . Since  X > NL,  it  follows  from  Theorem  2.13 
that  (2.I9)  may  be  rewritten  as 

+ P(0e(N)  + g(X)  t Z*-^  Cq(X,N,  s)  ) . (2.21) 

Also,  C^(X,N,s)  = zCQ(X,N,s),  so  the  first  equation  of  (2.16)  can  be 
rewritten  as 
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C°(X,N,s)  (1-^)'1Fq  + t^(l-OZ)-1  c\(K,K,s)  . 

Let  X X,  N N and  s = s > 1 be  fixed.  Denote  = C^fX^NjS). 
Then  (2.16)  - (2.18)  and  (2.21)  yield  a square  (N(s+L)  + 1)  nonsingular 
linear  system  that  can  be  solved  to  yield  the  costs  (C^)  for  the  fixed 

A A A ^ O 

values  X,  N and  s.  By  Corollary  2.14,  it  suffices  to  compute  only  C . 

Obviously  it  is  impractical  to  solve  such  systems  for  all  possible 

A /X  /v 

values  (X,  N,  s) , so  some  sort  of  parameterization  is  necessary.  Suppose 

AAA 

that  the  boundary  equation  (2 .2 L)  for  the  fixed  values  (X,  N,  s)  is 
replaced  by  an  equation 


N * N xv 

C-NL  -NL  + ^-NLfl  + ^Cc 


where  is  a dummy  variable.  Then  the  resulting  linear  system  is 

a a A 

independent  of  the  value  X and  solving  for  in  terms  of  will 

A A. 

A As  AA  A.  A X-NL  ' O j . S' 

yield  CQ  = a 4-  bC  ^ . Making  the  substitution  = g(X)  + Z C0  + 0e(*O 

/ 'Q  A A-  A.  A A V Ml  - 1 

then  yields  = ( a + bg(X) +bp^( N) ) ( 1 - bZ  ) . Since  this  is  true 

for  any  value  X chosen,  it  follows  that  C°(X,N,s)  = ( a +-  bgfX)  + bPfi(N)) 
xv  X-NL  - l 

( 1 - bZ  ) , for  all  X > NL.  Hence  using  the  dummy  substitution, 

A A 

the  system  need  only  be  solved  once  for  given  values  (N,s)  to  obtain 
the  functional  C^(',N,s),  which  can  then  be  minimized  over  [NL,  a>) . 

With  regard  to  the  parameter  N,  it  should  be  obvious  that  the 

/A  A 

respective  linear  systems  for  N and  N+l  would  be  very  similar.  In 

A A A 

fact,  the  first  N(s+L)  equations  of  the  system  for  N are  part  of 

A 

the  system  for  Nfl.  Hence,  it  seems  reasonable  to  expect  that  solving 
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for  the  functional  Cq(*,N,s)  should  aid  in  solving  for  the  functional 
C^(*,  N+ 1,  s) . The  next  two  sections  develop  procedures  for  recursively 
obtaining  the  functionals  g”(-,N,s)  far  N = 0,  1,  2,  ...  . 

With  regard  to  the  parameter  s,  recall  that  the  optimal  parameter 
value  s*  is  independent  of  X,  N and  g( • ) (Lemma  2.3).  Given  the 
recursions  of  the  next  section,  this  fact  is  used  to  develop  a 
simple  procedure  for  determining  s*  at  the  onset  of  the  computations. 
Once  s*  is  determined,  only  the  functionals  (C^( • ,N, s*) } need  be 
computed. 

The  remainder  of  this  section  uses  probabilistic  arguments  to 
investigate  the  functional  form  of  Cq(’,N,s)  that  can  be  anticipated. 
Let  C°(N,s)  denote  the  expected  discounted  costs  until  the  first  expan- 
sion, starting  from  spares  level  k > 0 (with  no  modules  initially 
connected),  when  N is  the  limit  on  temporary  module  usage  and  the 
constant  s-policy  with  parameter  s is  in  use. 

Lemma  2.1b. 

G°(N,s)  = Zk  C°(N,s)  , k>0  . (2.22) 

Proof.  Starting  from  spares  level  k > 0,  no  costs  are  incurred  until 
the  spares  level  first  becomes  nonpositive;  the  time  of  this  event  is 

-rTki  k 

given  by  and  E[e  ] = Z by  Lemma  2.15.  The  expected  costs 

beyond  T^,  until  the  first  expansion,  are  given  by  Gq(N,s)  (discounted 

relative  to  T^) . Hence,  G°(N,s)  = E[G°(N,s)  e = G°(N,s)  E[e  ^ J 

„k 


o, 


□ 


When  an  expansion  occurs,  all  temporary  modules  are  disconnected 
and  the  spares  level  changes  to  X-NL  > 0.  Hence  the  expected  inter- 
expansion costs  are  given  by 


X-NL 


t . | ^,X-NL  -0/M  . 

N,s;  2 c0(N’s) 


Theorem  2.  I7 . 


ft(N.s)  » g(X)  NL 


X > NL 


where 


0(N,s) 


z-(NLfl) 


> 0 . 


(2.23) 


(2.23a) 


Proof.  Starting  from  spares  level  zero,  the  time  until  the  first  expan- 
sion is  given  by  and  E[e  rfNL+l]  _ zNL+ 1 , by  Lemma  2.15.  At 

time  p an  expansion  of  size  X-t-1  occurs  at  cost  g(X)  and  the 

expected  costs  thereafter  are  c‘J  (X,N,s)  ZX'NL  C°(X.N,s)  (discounted 

X-NL  ’ ’ ' 0 ’ , 


relative  to  T, 


NL+1; 


Hence . 


C°(X,N,s)  = Gq(  N,  s)  + E((g(X)  + ZX‘NL  C°(  X,N,  s)  ) 

= C°(N,s)  + (g(X)  + ZX*NL  C°(X,N,s))  E[e'rlNL+l] 
= G°(N,s)  + (g(X)  + ZX-NL  C°(X,N,s))  ZNL+l 
= cJ(N,s)  + g(X)  ZNU1  + ZX+l  Cq(  X,N,  s) 
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NL  t-  L 


C0(X,N,s) 


U0(N,s)  + g(X)  Z 

T-1F1 

Gq(N,s)  z'(NLfl)  + g( X) 


,NL 


z-1  - ? 


r'.N.s)  , a(x)  znl 

i • r 


Z'L  - zx 


where  p(N,s)  is  given  by  (2.23a).  by  Assumption  2.2, 
and  d > 0.  Hence,  all  costs  incurred  are  nonnogative. 
and  it  then  follows  that  p(N,s)  > 0,  since  Z > 0. 


il  > TjL,  c > 0 
Thus,  G°(N, s)  > 0 

O 


Rearranging  (2.23a)  gives  c”(N,s)  P'N,s)  ZNI  ’ PiN,s)E[e 
Similarly,  using  (2.22)  gives 


f>(N,s) 


zNL+k+l 


E[e 


“rTNLfk+l  . 


k > 0 . 


When  the  initial  spares  level  is  k > 0,  the  time  of  the  first  expansion 
is  T._  . Tlius.  in  an  expectation  sense.  O(N.s)  is  the  equivalent 

lump  sum  payment  that  would  be  incurred  at  the  first  permanent  expansion 
time  in  lieu  of  the  incremented  charges  over  the  time  span  previous  to 
the  expansion.  Since  the  expected  inter-expansion  costs  are  G^  ^(N,s) 

= 0(N,s)  E[e  A+i]  and  the  inter-expansion  time  intervals  are  (in 
distribution)  of  size  p a similar  lump  sum  interpretation  of  P(N,s) 

can  be  given  with  regard  to  the  inter-expansion  costs. 
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Lemma  2.18.  p( • f s)  is  nondec reas ing,  s > 1. 


Proof.  The  expected  costs  over  the  time  interval  [0,  are 

G°(N,s)  " P„(N)  Z^'  l,  regardless  of  whether  the  module  limit  is  N 
or  1.  The  expected  costs  over  the  greater  (with  probability  1) 
time  interval  [0,  T(N+1jl+1)  are  G°(N+l,s)  - 0e(N+l)  Z!Nfl)Lfl. 
Hence,  using  (2.23a), 


since  the  disposal  costs  Pe( ‘ ) are  assumed  nondecreasing.  □ 

By  Theorem  2.17,  the  coefficient  p(N,s)  completely  determines 
the  functional  C°(*,N,s).  Hence,  the  next  section  focuses  on  developing 
a procedure  for  recursively  computing  p(N,s)  over  N = 0,  12  ...  . 
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2.7.  Recursive  Computation  ol'  the  Functionals  { • , N , s) , N > 0 } 

As  demonstrated  in  the  previous  section,  the  coefficient 
0(N,s)  completely  determines  the  functional  Cq(«,N,s).  For  a given 
value  s > 1 this  section  provides  a procedure  for  recursively  com- 
puting the  coefficients  {p(N,s),  N 0,  1 ...)  and  hence,  the 

functionals  {C°f*,N,s),  N - 0,  1,  ...}.  The  next  section  provides  a 
procedure  for  determining  the  optimal  disconnect  parameter  s*  at 
the  onset  of  computations.  Thus,  the  algorithms  presented  here  need 
only  be  implemented  for  the  given  value  s = s*. 

Definition  2.%  For  m > 0,  let  A(m)  denote  the  square  tridiagonal 


matrix  of  dimension 

m+1 

with  nonzero  elements  given 

by 

> 

H* 

II 

1 , 

i = 0,  1,  • • • . 

m 

A . . . ( m)  = 
1,1+1  ' 

i = 0,  1,  ..., 

m-  1 

A.  . ,(m)  = 
l.i-l  ' 

-a  , 

i = 1,  2,  ..., 

m 

Definition  2.6.  Let  C,  R C TRmf  ^ and  Q t IR1  R*  (m+2).  The  square 
linear  system 


is  a Q- system  if  Q is  nonsingular  and 
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1 u 

V A 


where  A = A( m)  f U1  eBnul)  V t !Rmfl, 


and 


t 


j / e 


j = t 


i = 0 


i { 0 


for  some  af  b and  some  £ > 0. 

The  following  lemma  demonstrates  that  the  solution  to  a 
is  determined  by  A~ ^ . 


Lemma  2 . I9 . If  the  square  linear  system 


is  a Q-system,  then 
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Q-system 


1 


J 


C-  = ( 1 - abA*  l)  ' 1 ' R - bA"  1 R) 

0 io'  0 C.  ' 


C - ac0  A-*  . 


Proof.  Partition  Q as 


v a 


where  e t H,  UT  ( F^1,  V < 1 and  Q f R(n^l)x(ra+1) . from 


Definition  2.6,  Q Q ^ I implies  that 


e + au=  1 


eu  + UA  - (0,  0) 


1 2.27 ) gives,  using  Definition  2 .6, 


U - -elJA'1  = -ebA",1  . 

4/  • 


Thus.  -ebA,Q.  Substituting  this  result  into  (2,26)  gives 


€ = (i  - abA^V1 
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Hence,  from  (2.28), 


U - -d  - .bA^1)-1 


and  thus, 


C 


0 


o - abA«>  <ro  • M«!  s> 


Using  Definition  2.6, 


aCQ  e(  0)  + AC  = R , 


where  e(0)  = (1,  0,  ...,  0)T  ]Rm+1.  Hence, 


C = A_1(R  - aCQ  e(  0) ) 


A' bR  - ac„  A'* 


The  next  lemma  treats  the  special  case  N 


Lemma  2.20. 


where 


= ^a(°,s)  ♦ Pc(0)  , 


Proo^.  When  N = 0,  the  transition  equations  aie  independent  of  s 
temporary  facilities  are  never  used) : 


C 


o 

k 


Zk  C° 

0 » 


k > 0 


-QC?  + C°  = F°  «•  0(g(X)  + pe(0)  + C°  ZX)  . 


Substituting  Cx  _ ZC°  in  the  second  expression  and  collecting  terms 
yields 


co(x,  0,  ■) 


P'1  1-0  + tjo)  + g(X) 

0_1(i-az)  - zx 


F0  + 0e(O)  } + 8fX) 


The  last  equality  follows  since  Z satisfies  aZ2  - Z +.  0 = o,  by 
Theorem  2.3.  Thus,  p(0,s)  *fl(0,s)  + PJo),  where  p^O.s)  & 1 F°, 

independent  of  s > 0. 


Denote 


F"  = 


(F 


-(n-l)L+s-l>  F-(n-l)L+s-2» 


F°  t)T 

- nL' 


IK 


s+L 


Let  (e(j))  denote  the  unit  vectors  (where  the  dimension  of  e(j)  is 
determined  by  the  context  in  which  it  is  used);  that  is,  (e(j)).  1 

and  (e(j))1  = 0,  i / J,  for  j = o,  1,  2,  ...  . 


( since 
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Theorem  2.21, 


and  £>a(0,s) 


Let  s > 1.  Denote  m = s+L-1 
p"1  F°.  Then,  for  N > 1, 


A A (in) 


of  0,  s)  = 0 


0(N,s)  0a(N,s)  + 0e(N)  , 


( 2.  29) 


where 


0 (N,s)  l—  10  ( N- 1 s)  +A'L(FN  fao(N-l,s)  e(  0) ) ) (2.30) 

a PA'1  3 S> 

sm 

and 


o(N,s)  = a:1  [FN  + Oo(N-l,s)  e(0) 


p0a(N,s)  e( m) ) 


(2.31) 


Proof.  The  proof  is  lengthy  and  is,  therefore,  only  outlined  here. 

See  Theorem  3.21  of  [9]  for  details. 

The  proof  was  the  "dummy  variable"  substitution  previously 

discussed  in  Section  2.6.  For  N = 1,  the  transition  equations  form 
a Q-system  of  dimension  m + 2.  Recall  (Section  2.6)  that  the  first 
N(s+L)  equations  of  the  linear  systems  for  N and  N + 1 are 
identical.  In  general,  the  proper  substitution  of  the  results  for  N 
into  the  transition  equations  for  N + 1 will  similarly  yield  a 
Q-system  of  dimension  m + 2.  The  theorem  then  follows  by  induction 
on  N using  Lemma  2.23. 
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Theorem  2.21  provides  the  necessary  recursion.  However  in  order 
to  implement  the  recursion,  A ' and  A^^  must  first  be  computed 


(where  A A(s+L-1).  The  next  lemma  provides  a simple  means  for  computing 
1. . ’ 


A'1  ,,  A'1 

A and  A. 


s . 


•1 


Lemma  2.22.  A (0)  = [Ij  and  for  m > 0, 


A"*(jir1)  = (A'\m)  + Q05(m)  A"Vm)  A‘L(m),  3&(m)  AT^(m)), 

i = 0,  1,  ...,  m (2.32) 

and 

Am!l  = (Q6(tn)  A‘ \m) , 6(m))  , (2.33) 

where 

&(m)  = (1  - QpA'j^m))'1  . (2.34) 


Proof. 


By 


Definition 


A(m+l)  = 


2.5,  A( 0)  = 1 1] 

A(m)  B 

C 1 


A (0);  alsOj 


where 


(2.35) 
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B = (0,  0,  ..., 

0,  -f3) 

T C B^1 

(2.35a) 

c = (0,  0, 

0,  -a) 

. CT  Rm+1  _ 

(2.35b) 

Partition  A \m+l)  as 

-i  r a 

A (m+1)  -- 

B ' 

> 

(2.36) 

_ (* 

&(m)_ 

, _ _.(m+l)x(m+l)  _ 

whore  Q . 3R  ' ' , B 

mma. 

• Bmfi,  and  6(m)  B. 

Since 

A(m+1)  A (m+1)  = I,  it  follows  from  (2.35)  and  (2.36)  that 

A(m)  G + B£ 

1 , 

(2.37a) 

A(m)  3 + 6(m)B 

= (o,  • 

o)T  , 

(2.37b) 

ca  + A 

= (o,  . 

o)  , 

(2.37c) 

C 3 + &(m) 

= 1 . 

(2  .37  d) 

From  ( 2 . 37b) } 

B = -b( m)  A" l(m)B  . (2.38) 

Substituting  (2.38)  into  (2.3?d)  gives 
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1 


-6(m)  CA  ^(m)B  + 6(m)  = 

* "j(  m)  - (1  - CA”  l(  m) B) " ‘ . (2.39) 

Using  (2.37  a), 

A * m)  B -pA  / . (2.40) 

Hence,  using  2.3  b) 

CA  *('n)B  ^ CtpA  ' ( ni) 
mni  ' 

Substituting  the  above  expression  into  (2.39)  gives 

b(m)  = (1  - at3A^(m))"1  , 

which  verifies  (2.34). 

Using  (2.37a), 

a A~L(m)  - . (2 .4 1) 

From  (2.37c),  (2.39)  and  (2.4l), 


jL 


b6 


CQ  t f C(A_l(m)  - A"  1 ( m)  BC)  + ft 

CA  *(m)  - ( CA  ^(m)B-l)l*  0 

CA_1(m)  = (CA*1(m)B  - l)!4 

ft  - (CA"l(«)B  - l)'1  CA"l(m) 

= -6(m)  CA  ^(m) 

Using  ( 2 . 35b) f 

CA_1(m)  = -CiA'l<m)  . 

J m.  1 

Substituting  (2.43)  into  (2.42)  yields 

ft  afc(m)  A-  l(  m)  . 

m , 

Substituting  (2.40)  into  !2.38j  similarly  gives 

6 - p&(w)  A_1(m)  . 

7 .nr  7 


(2.42) 


(2.43) 


(2.44) 


(2.45) 


Using  (2.40)  and  (2.44)  yields 
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A V m)  B<^  ■ tiA  V m)  6(m)u  A ^(  m) 

. m m. 


-UJ3b(m)  A Vn)  A Vb) 
' . m ' m.  ' 


Substituting  the  above  expression  into  (2.41)  gives 


Q A ^ ; ra)  t Gp  (m)  A ^(m)  A^(tn) 


(2.46) 


Using  (2.  16,  ( (2.32)  follows  from  (2,45)  and  (2.46);  (2,31) 
follows  similarly  from  (2.44). 


Using  Lenina  2.22,  a simple  algorithm  for  computing  A,"  (st-L-1) 

.1  ’ ’ 
and  A^  (s-tL-l)  can  be  written.  Let  m^  min(s,  L-l)  and  nr,  = max{s,L-l), 

Then  the  algorithm  can  be  outlined  as  follows: 

(a)  Compute  (m^)  using  (2.33),  for  m = 0,  ...,  m^- 1. 

(b)  Compute.  A 1 (m,,)  using  (2.32)  and  A'1  (m0)  using  (2.33),  for 

***  | • 1 HI,  . a 4l 

•*  j nVj"  1 • 

(c)  Compute  A"1  ( s+L- 1)  and  A-1  (st-L-1)  using  (2.32)  for 

m\9  mo* 

m = nu,  . . . , s+L-2. 


It  should  be  noted  that  (2.33)  must  also  be  used  in  part  (c)  in  order 

to  maintain  A 1 ( m)  for  in  m, SfL-2.  The  actual  algorithm  is 

m.  . 7 7 ° 

given  below  ind  uses  three  vectors  B1  , i 1,  2,  J,  of  arbitrary  dimension. 


r 


1 -I  f''  — 1 

Algori thm  B^  (Compute  B A (s+L-1)  and  B Aj^  j (s+L-l)). 

( 1)  If  s >L-1,  go  to  step  (2);  otherwise,  i <-  1,  j «- 2,  m^  ♦-  s, 
«-  L-l,  m «-  0,  m.  s-tL- 1,  B 1 • ( 1)  and  go  to  step  (3)  . 


m. 


2 

i «-  2. 


j «-  1,  m j «-  L-l,  «-  L-l,  m <-  0,  m^  *-  Si-L-  1,  B «-  ( I)  . 


(2) 

(3)  If  m >1^-1,  go  to  step  (5);  otherwise,  continue. 

(4)  & *-  ( 1 - OpB^)  , B*  (Qf  B1,  t) , me-m+1;  go  to  step  (3) 

(5)  BJ  *-  B 1 . 

(6)  If  m >111,-1,  go  to  step  (8);  otherwise,  continue. 


(7)  6 - (1  - Q^B j ) ' 1 , 


m «-m+l;  go  to  step  (6). 
(8)  B:  <-  Bj. 


(B1  + OpbB^  P6B^),  B ' *-(a£>Bj,  &), 


(9)  If  m >m^-l,  stop;  otherwise,  continue. 
(10)  6 «-(  1 - OfJBj'1.  B 


nr 


(B1  +■  Q66B1  B^.  (36B1)  , B^  <- ( B"  + QpbB"  , 

1 m ’ nr  * m ’ 

p 5 < 

PbB^) , B *- (dfcB  , &),  m t-m+1;  go  to  step  (9). 

When  the  algorithm  terminates  in  step  (9),  vector  B^  is 
1 2-1 

A"  (s+L-1)  and  vector  B is  A (s+L-1).  Given  these  vectors,  the 
S . L- 1^  , 

recursion  of  Theorem  2.21  can  then  be  used  to  compute  the  coefficients 
0(  - ,s). 


Algorithm  B,  (Compute  0(N,s),  N = 0,  1,  2,  . . . , N) . 

12 

( 1)  Compute  vectors  B and  B~  (use  Algorithm  B^) . 

(2)  0a  <-  P'  ' F°,  a 0,  m «-  s+L-  1,  N <-  0,  a <-  ( PB*)  1 . 

(3)  ^(N,s)e()a<Pe(B). 

(4)  If  N > N,  stop;  otherwise,  continue. 

(5)  N *-  N+l,  B5  e-F**,  Bo  -Bq  + - a(*>a  + b1'B'>»  Bm  «“  Bm  “ 

2 5 

a ■&;  go  to  step  (3), 
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The  vector  B^  used  in  step  (5)  is  a work  vector.  The  calculations  in 

step  (5)  follow  the  recursions  (2.30)  and  (2.31),  given  B1  = A*  *"(  s+L- 1) 

s. 

2 - 1 

and  B = ^ ^ (s+L-1)  from  Algorithm  B^. 

Algorithms  B^  and  provide  a means  for  determining  the  functionals 
(C°(-,N,s),  N > 0)  for  a fixed  given  value  s > 1.  By  Theorem  2.11, 
only  the  functionals  (C°( • ,N, s*) , N > 0)  need  be.  computed.  Therefore, 
the  next  section  focuses  on  the  problem  of  determining  the  optimal  dis- 
connection parameter  s*  prior  to  initiating  Algorithms  B^  and  B^. 


2.8.  Determining  s* 

In  this  section,  a procedure  is  derived  for  computing  s*,  the 
optimal  disconnection  parameter.  By  Lemma  2.3,  the  value  s*  is  independent 
of  X,  N,  Tj  and  g.  Therefore,  let  X = L,  N = 1,  = 0,  g = 0 and 

denote 

Cq(8)  = c°(l,i,s;  8 5 °,  rl  = 0),  S >0 

The  only  costs  contributing  to  C°( s)  are  the  charges  c,  d,  it  and 
0g(  1) . f)  ( 1)  is  incurred  only  at  expansion  epochs  and  those  epochs  are 
independent  of  the  disconnect  parameter  s.  Therefore,  only  the  temporary 
module  costs  resulting  from  connection  (c),  disconnection  ( d) , and  rental 
(it)  charges  vary  with  s.  That  is,  the  difference  E°(s+1)  - C°(  s) 
must  be  a positive  affine  multiple  of  -H(s+l,s),  s > 0.  Hence  it 
follows  from  Definition  2.h  that 
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(2.47) 


t 


I 


Cq(s*)  = min  C°(s)  . 

s > 0 

Fur thermore^  it  foilows  from  Theorem  2.8  that  ~^q(')  is  strictly 
unimodal  about  [s*-l,  s*j. 

The  results  of  the  last  section  provide  expressions  for  C°(s)? 
s > 1.  When  s = 0,  Model  II  reduces  to  Model  I.  The  next  lemma 
provides  a means  for  identifying  instances  where  s*  = 0;  in  these  cases; 
Model  I will  be  applicable.  Proceeding  as  before^  the  transition  equations 
for  the  case  s = 0 and  N = 1 are 


cJ(X,l,0)  - e(l-oz)-1  C^X^O)  = ( l-QZ)"1  F°  , 

-QcJ+i(x,I,°)  + CJ(X,1,0)  - X(X,  1,0)  = F*,  -L+l  < k < -1 
-ccl_M(x,i,o)  + c_Ll(x,i,o)  . f_1l  + p(K)e(i)  + g(X)  + ZX‘L  CJ(X,  1,0))  . 


Lemma  2.23.  Denote  A = A(L-l).  Then 


Cq(X,1,0) 


0,(1, o)  ♦ 0JO  + g(x) 


z-1  - zx 


where 


0ad,o) 


P*. 


-1 
0,L-1 


(P 


-1 


v1 


F1) 


Proof . See  Lemma  3.23  of  | 9 ] . 


(2.48) 


( 2.49) 
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Denote  P ( s)  f\  ( 1 s;  if.  - 0)  s > 0.  Then  by  Theorem  2.21 

cl  cl  J.  * 

and  Lemma  2.23, 


pac  *)  +V1) 

z-1  - zL 


s > 0 . 


( 2 . 50) 


Hence,  by  (2.4?) 


*>(■*)  = min  PS*)  , (2.51) 

s > 0 

and  furthermore  -f>  (•)  is  unimodal  about  [a*-l  s*].  Using  Theorem  2.21 

and  Lemma  2.23,  a recursive  procedure  will  be  derived  for  evaluating 
0 (s),  s = 0,  1,  ...}  in  order  to  determine  s*. 


Lemma  2.  24.  A"*(0)  = [lj 


and  for  m > 0, 


A'^m+i)  = (b'(m),  06'(m)  A’^m))  , (2.52) 

and 

A‘Vm+1)  = (Qfc'(m)  A^^m),  A*1^  (m)  + 0p6'(m)  A"^  Q(m)  A^l(tn)), 

i = 1,  . . . , m+1  ( 2.  53) 

where 

5’(m)  = (1  - apA^(m))-1  . (2.54) 


^roof.  By  Definition  2.5f  A(0)  = flj  A_l(0);  also. 


A(m+  1) 


1 B 
C A(  m) 


(2.55) 


where  B = (-0,  0,  ...,  0);  BT  6 *mfl  .and  c = (_a>  0>  ^ 0)T  £ Rn 

Partition  A (nw-1)  in  a manner  analogous  to  (2.55).  The  lemma  then 
follows  as  in  the  proof  of  Lemma  2.22  using  A-1(m+l)  A(m+1)  « I.  See 
Lemma  3. 24  of  [9]  for  details. 


Definition  2.Y.  For  s > 0,  let 

aL(s)  = A^(.+L-l)  , 
a2(s)  = A^(s+L-l)  , 

a3(s)  = Ad,1s+L-l(s+L"l>  » 

a4(s>  “ » 

ei(s)  = A8^S+L-L)  i > 

%(*)  = A”  X(  a+L- 1)  l , 

who  re 

1 = (1,  1,  1)T  • 


Lemma  2.25.  For  m > 0, 


and 


\i(»D  • »■(•)  , 

‘o'nwlf”*1)  **'<■>  *^(») 
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(2-56) 


(2.57) 


(2.58) 


Aq‘  ( nif  1)  = h'  ( m)  (1  + (m)^}  , 

where  &' (m)  is  given  by  (2.54). 

Proof.  The  results  follow  directly  from  (2.52)  and  (2.54)  of  the 
previous  lemma.  □ 

Lemma  2,26.  For  s > 0f 

a^s+1)  * (1  - opa^s))'1  , 
a2(s+l)  = aa^s+1)  a2(s)  , 
a^(  s+1)  = fia^s+1)  a^(  s)  , 

(2.59) 

ai;(  s+l)  = a[f.( s)  + 3a2(s-t-l)  a^(s) 
e1(s+l)  = e^s)  + a2(s+l)  {1  Pfi^s)} 
e^s+l)  = a1(s+l)  {1  + 36^  s)  J . 

Proof.  Using  Definition  2.7  with  Lemmas  2.24  and  2.25, 

a^a+l)  = A^(s+L)  = &"(s+L-l)  = (1  - QpA^(  s+L- 1)  )* 1 

= (1  - Opa^s))'1 

a2(8+l)  = a;Jj  q(  s-t-L)  = Q&'(s+L-1)  A‘^(s+L-1) 

- aa^s+l)  a2(s) 
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s + 1) 


\(  s+l) 


eL(  s+i) 


Cgfs+i) 


The 


* A0,‘s.L<  s*1)  S8'(».L-1) 

= pa1(s+ 1)  a (s) 

' As!l,stL<s+L> 

* ‘s^L-l'8*1-1)  * ( s+L-  1)  A*  s+L-1) 

= au(s)  + a2(s+l)  ?AQ1g+L_ i( S+L~ !) 

= \(s)  + ^a2(s+l)  a^(s) 

- A's!i,.(stL>  i 

= 0&’(  s+L-1)  A~sq  (s+L-1)  + A”  X(  s+L-  1)  l 

+ QpS' (s+L-1)  A~sq  (s+L-1)  A‘L(  s+L-  1)1 
= eL(s)  + 06*  (s+L-1)  a;  J(  S+L-1)  (1  + peg(s)} 

= 61(s)  + a2(s+l)  U + (^(s)) 

= A~  \ s+L) 1 

= &' (s+L-1)  + p&'(  s+L-1)  A0^(  s+L-  1)1 
= 6’ (s+L-1)  {1  + s+L- 1)1) 

= ^(s+l)  {1  + P£2(s)  ) . C 

significance  of  Definition  2.7  is  given  by  the  following 


theorem. 
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Theorem  2.2  7. 


0a(s)  = c + (pa4(s))'  {c  r Od  a,,(s)  + iiei( ») } , s >0  . (2.60) 


Proof.  By  Theorem  2.21  and  Lemma  2.23, 


5*(s)  ‘ T1 — L~ ~r  ls'1  Fo * As.1(,+L-1)  pl) 


(2.61) 


Given  that  = 0, 


F0  = 


1 c , T 

f = (od  + ii,  n,  ii,  ...,  ii,  tic  + ii)  eiR 


Substituting  these  expressions  into  (2,61)  and  Definition  2.7  yields 


8a<“>  - 5T7tT  ‘c  * *";Wl)l) 

* * 9 


'^a'Vs7  ic  + ada2^s)  h tJca4(s)  +n-e1(s)) 


c + (pa^s))’1  {c  + Qda2(s)  + 11^(8)) 


The  desired  recursion  can  now  be  outlined  as  follows: 
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(a)  Use  (2.56)  - (2.58)  to  compute 


a^fO)  = a2(0)  = A00(L-1)  , 
a3(0)  = a^O)  = A-|L1(L-1)  , 

and 

e^O)  = £2(0)  A-^(L-l)l  . 

( b)  Recursively  compute  p^(s),  s > 0 using  (2.60)  with  (2.59)  until 
s*  is  identified  by  0 ( s*)  < (s*+l). 

£L  SL 

The  relationship  identifying  s*  in  part  (b)  follows  from  the  unimodal 
property  of  -f>  (').  The  algorithm  is  stated  below. 

Algorithm  B/(  (Find  s*). 

(1)  aL  1,  a5  «-  1,  <-  1,  n>  «- 0. 

(2)  If  m > L- 1,  go  to  step  (4);  otherwise,  continue. 

(})  aL  4-  ( 1 - Cipaj)'1,  a^  4-  pa^,  *-3^1  + pe t),  m <-  m+1;  go  to 
step  (2). 

(i+)  s*  4-0,  a2  <-  a,  a^  *-  ^ 4-^,  pfl  *-(Pa[))'1  (c  + oda2  + n^). 

(5)  aj_  e- (1-Qpa^-1,  a2  ♦- Oa^,  a^  a4  + pa^,  a^  «-  v «-  1+PEg, 

£1  <_gl+  a2v>  £2  *“alv>  ^a  (Pa',)'l(c  + ^2  + ,rei)- 

(6)  if  stop;  otherwise  p 4- $ s*  4- s*  + l,  and  go  to  step  (5). 

Steps  ( 1)  - (5)  of  the  algorithm  compute  a^(0),  a^(0)  and  £^(0)  using 
(2.56)  - (2.58).  Step  (1<)  computes  $ (0)  - c.  Step  (5)  successively 
computes  $a(  s)  - c,  s > 1,  and  step  (6)  compares  ^a( s)  - c with 
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$a(s"l)  - c.  The  algorithm  terminates  since  s*  < +oo^  by  Leimna  2.2. 
When  termination  occurs,  the  optimal  disconnection  parameter  is  given 
by  s* . 


2.9.  Summary  and  Statement  of  the  Expansion  Size  Optimization  Problem 
Let  k = {0,  L,  2L,  . ..,  NL},  where  N denotes  an  upper  bound 
on  temporary  module  usage  based  on  physical  considerations.  For  K £ k, 
denote 


0(K)  = P( K/L,  s*)  and  Cfc(XfK)  = C°(X,  K/L,  a*),  k >0  . (2.62) 

Given  K t k and  an  initial  spares  level  kQ  > 0,  an  optimal 
permanent  expansion  size  X*(K)+1  for  Model  II  is  given  by 


Cq(X*(K.),  K)  = min{C0(X,K)  : integer  X > K)  , 


where 


VX,K) 


- m ± &{*)  zk 

z-1  - zx 


X > K . 


£(•)  is  nonnegative  and  nondecreasing  over  k; 

0 < Z = — 20^^”  " given  a and  (3 

satisfying  (2.15)  (QZ2  - Z + f}  = 0)  ; 


(2.63) 


(2.64) 


(2.65) 


(2.66) 
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r 


I 

Cr(X,K)  ZkCQ(X,K),  k > 0 ; (2.67) 

and 

— t‘T 

zk  =,  E[e  kJ,  k > 1 . (2.68) 

(2.64)  follows  from  Theorem  2.17,  (2.65)  follows  from  Theorem  2,17  and 

Lemma  2. 18,  (2.66)  and  (2.67)  follow  from  Theorem  2 . 13 ^ and  (2.68) 

follows  from  Lemma  ^.15.  The  fact'  that  X*(K)  is  given  by  (2.63), 
independent  of  k^  _>  0,  follows  from  Corollary  2. 14. 

Algorithm  B. , presented  in  the  previous  section,  provides  a 
simple  means  for  determining  s*.  If  s*  = 0,  then  Model  I of  [7] 
is  applicable.  If  s*  ^ 0,  then  Algorithms  and  B(J  of  Section  2.7 
provide  recursions  for  computing  the  coefficients  0(K)  = 0(K/L,  s*) 
over  K c k.  Once  these  coefficients  are  known,  the  task  of  determining 
X*(.)  over  k becomes  the  series  of  single-var iable  minimization  problems 
given  by  (2.63)  under  conditions  (2.64)  - (2.68).  Notice  that  this 
minimization  problem  is  identical  to  that  resulting  for  Model  I of 
[7).  The  optimization  problem  is  treated  in  [8J.  By 
Corollary  2. 14,  the  optimal  temporary  facilities  usage  limit  K* 
(alternatively,  the  optimal  temporary  modules  limit  N*  = K*/L)  and 
the  associated  optimal  expansion  size  X*+l  are  then  given  by 

C0(X*,  K*)  = min{C0(X*(K),  K)  : K £ k)  . 


jL 
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It  is  shown  that  an  iptimal  module  removal  policy  has  the  form:*  ^remove  a 
model  only  if  at  least  s*  additional  units  of  unused  capacity  will  remain 
available''. * A simple  algorithm  is  given  for  determining  the  optimal 
parameter  s*.  jk 

For  a givenXlimit  (K)  on  temporary  facility  usage,  the  form  of  the 
expected  discounted  cost  functional,  parameterized  in  the  expansion  size 
(X+l) , is  derived.  Recursions  are  given  for  determining  these  functionals 
over  all  feasible  values  for  K.  A sequel  paper  [8]  treats  the  problem  of 
minimizing  the  functionals  in  order  to  find  optimal  expansion  sizes. 


